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ABSTRACT 


Fourier  Trans form  Spectroscopy  has  been  used  to  study  the 
infrared  emission  from  a premixed,  methane  - air  flame  using  a 
flat  flame  burner,  in  the  region  1700cm  1 to  7000cm-^.  Some 
preliminary  data  was  also  obtained  in  the  region  from  400cm  1 
to  1700cm-1 • All  spectra  were  taken  with  1cm-1  resolution  or 
better.  Spectra  were  obtained  at  various  heights  above  the 
burner  base  and  under  a variety  of  mix  ratios.  From  the  rota- 
tional structure  in  the  CO  fundamental  band,  rotational  and  vi- 
brational temperatures  were  obtained  for  a rich  fuel  mix. 
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ANALYSIS  OF  COMBUSTION  PLASMAS 
BY  FOURIER  SPECTROSCOPY 


I.  Combustion  Analysis 

A.  Objective 

The  purpose  of  this  research  is  to  establish  the  fea- 
sibility of  using  Fourier  Transform  Spectroscopy  to  analyze 
combustor  plasmas  of  AF  interest.  The  need  to  advance  the 
quantitative  knowledge  of  combustion  processes  stems  from 
a combination  of  reasons.  There  exists  a finite  and  exhaus- 
tible supply  of  conventional  fossil  fuel,  dictating  the 
requirements  to  maximize  the  efficiency  of  converting  this 
fuel,  as  well  as  to  develop  alternate  fossil  fuels  such 
as  derivatives  of  shale  oil  and  coal.  Coupled  with  these 
requirements  are  existing  and  projected  Federal  regulations 
regarding  exhaust  emission  in  the  troposphere  and  strato- 
sphere. 

B.  Background 

Theoretical  modeling  of  the  combustion  plasma  relies 
on  a knowledge  of  the  species  in  the  plasma  and  their 
respective  reaction  rates.  Basic  data  of  this  nature  is 
still  inadequate  however.  Tools  for  combustion  diagnostics 
include  probes,  such  as  mass  spectrometers,  and  several 
optical  techniques.  Mass  spectrometer  measurements  yield 
data  on  the  identity  of  the  species  and  their  relative  con- 
centration, but  the  technique  is  limited  by  its  ability 
to  sample  only  reaction  products  downstream  from  the  source. 
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and  by  the  fact  that  a variety  of  hydrocarbons  have  the 
same  mass  number.  Two  powerful  optical  techniques  have 
recently  been  developed,  but  similarly  do  not  yield  a com- 
plete description  of  the  physics  of  combustion.  The  laser 
raman  scattering  (LARS)  signal  from  the  combustion  chamber 
is  used  to  obtain  spatially  resolved  profiles  of  the  major 
species  and  temperature  of  the  combustion  plasma.  This 
technique  has  limitations  due  to  a lack  of  signal  strength. 
Coherent  antistokes  Raman  spectroscopy  (CARS)  offers  a 
significant  advantage  over  incoherent  Raman  spectrosocpy 
at  atmospheric  pressures  in  that  the  signal  strength  are 
several  orders  of  magnitude  larger  than  spontaneous  Raman 
scattering  for  the  same  target  density.  While  both  the 
CARS  and  LARS  requirements  for  high  peak  power  tunable 
laser  have  been  met  in  part  by  advances  with  tunable  dye 
lasers,  the  applications  are  still  limited  to  cooperative 
situations. 

The  emission  and  absorption  of  radiation  associated 
with  the  combustion  plasma  contain  detailed  information 
about  the  combustion  species  and  their  reactions.  Since 
the  frequency  of  the  radiation  emitted  or  absorbed  by  a 
molecular  gas  is  dependent  on  the  energy  level  spacing  of 
the  molecule,  and  the  relative  intensities  of  the  radiation 
is  a function  of  the  energy  level  population  distribution. 
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the  identity  of  the  emitting  or  absorbing  gas  and  its 
temperature  can  be  obtained  from  analysis  of  the  emission 
or  absorption  spectrum.  The  intensity  of  radiation  in 
the  infrared  portion  of  the  spectrum  resulting  from  vibra- 
tional and  rotational  transitions  is  very  weak,  however, 
so  that  emission  and  absorption  analysis  have  been  applied 
most  successfully  in  a laboratory  environment.  Now,  with 
the  implementation  of  Fourier  transform  spectroscopy  (FTS) , 
the  sensitivity  of  infrared  spectrometers  have  been  increased 
by  up  to  100  times  the  sensitivity  of  dispersive  systems. 

The  power  of  FTS  is  that  it  is  much  faster  than  dispersive 
systems,  and  this  speed  may  oe  used  to  improve  the  signal- 
to-noise  ratio  by  coherent  addition  of  signal.  In  addition, 
recent  technical  developments  with  data  handling  systems 
have  greately  enhanced  the  capabilities  of  FTS,  allowing 
for  rapid  signal  averaging  and  conversion  from  time  to 

^frequency  domain,  as  well  as  sophisticated  manipulation 

t 

of  data  stored  in  the  system.  FTS  has  been  successfully 
applied  to  studies  of  infrared  spectra  of  aqueous  solutions, 
kinetics  of  chemical  reactions,  and  infrared  chemiluminescence. 
Applied  to  combustion  processes  the  FTS  system  assembled 
here  has  the  capability  to  identify  transient  and  product 
species  in  the  combustion  plasma,  as  well  as  the  relative 


concentration  and  temperature  of  the  species.  This  will 
be  demonstrated  in  the  results  for  rotational  temperature 
of  a flue  flame  burner. 

II.  Fourier  Transform  Spectroscopy 

A.  Introduction 

Fourier  Transform  Spectroscopy  (FTS)  is  a multiplex 
technique  for  measuring  the  infrared  power  spectrum  of  an 
emitting  source.  Multiplexing  implies  than  many  separate 
streams  of  information  are  simultaneously  transmitted  along 
a single  channel,  and  in  FTS,  information  about  a large 
number  of  spectral  elements  is  simultaneously  gathered 
through  a single  detection  channel.^-  In  an  ordinary  spectro- 
graph, information  about  many  spectral  elements  is  gathered 
simultaneously  by  the  photographic  emulsion;  however,  each 
emulsion  grain  acts  as  a detection  channel.  Thus  the 
ordinary  spectrograph  is  not  a multiplex  system,  it  is  more 
properly  a multichannel  system.  In  multiplex  spectroscopy 
the  intensity  of  each  spectral  element  of  a broad  spectrum 
is  coded,  energy  from  the  entire  spectrum  detected,  and 
then  the  spectrum  decoded.  The  coding  used  in  FTS  is  the 
orthogonal  functions  sine  and  cosine.  Specifically,  the 
intensity  of  each  spectral  element  6a  (6a  in  cm  ^ is  the 
resolution  of  the  instrument)  is  sinusoidally  modulated 
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at  a frequency  which  is  proportional  to  the  wavenumber 
of  the  element.2  An  advantage  to  this  approach  lies  in  the 
fact  that  since  information  about  the  entire  spectral  region 
is  gathered  simultaneously,  there  is  a savings  in  measure- 
ment time  over  a sequential  measurement  technique,  such 
as  a normal  monochrometer.  To  see  this,  let  the  spectral 
region  of  interest  extend  from  to  a 2»  Then  the  number 
fo  spectral  elements  in  the  region  is  given  by; 

°2~°1 
m = ~ 

6a 

and  the  time  to  sequentially  measure  the  spectrum  is  mT, 
where  T is  the  time  devoted  to  the  measurement  of  a single 
element.  In  the  multiplex  technique,  the  time  required  to 
measure  the  entire  spectral  region  is  just  T.  If  the  measure- 
ment is  detector-noise  limited  (the  detector  is  the  chief 
course  of  noise)  then  the  signal-to-noise  ratio  associated 
with  the  measurement  of  a single  element  goes  as  the  square 
root  of  the  measurement  time  for  the  element.  If  one  requires 
that  the  total  measurement  times  for  the  two  techniques 
be  equal,  say  T,  then  the  time  devoted  to  each  element  in 
the  sequential  technique  is  T/m,  while  the  time  devoted  to 

each  element  in  the  multiplex  technique  is  T.  In  the 

1/2 

first  case,  the  signal-to-noise  ratio  goes  as  (T/m) 
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Thus  for 


while  in  the  second  case  it  goes  as  (T)^/^. 
a fixed  total  measurement  time,  the  multiplex  technique 
results  in  a signal-to-noise  improvement  of  over  a 

sequential  measurement  technique.  This  reduced  measurement 
time  or  signal-to-noise  improvement  is  known  as  the  multiplex 
of  Fellgett  advantage . ^ ' 2 ' ^ ^ 

B.  Basic  Equations 

To  realize  the  multiplex  advantage,  the  spectral 
information  from  the  source  must  be  coded  and  detected. 

The  device  that  accomplishes  this  in  FTS  is  the  Michelson 
Interferometer  (or  one  of  its  variants) . Figure  II 
illustrates  an  interferometer  and  its  collimating  optics. 
Radiation  incident  on  the  beam  splitter-compensator 
is  amplitude  divided  into  two  parts.  These  disturbances 
travel  to  mirrors  and  M2  and  are  recombined  by  the  beam- 
splitter-compensator after  traveling  over  two  different 
paths  X^and  X2.  The  resultant  disturbance  is  then  focused 
on  the  detector.  One  of  the  mirrors  is  movable  so  that 
the  path  difference,  X2~X1=X,  can  be  varied.  Consider  the 
intensity  that  reaches  the  detector  at  a fixed  path  difference. 
Before  doing  so,  consider  the  source  to  have  the  following 


BASIC  MICHELSON  INTERFEROMETER 


M 


I 


Figure  II-I 


(2) 


it  is  ergodic,  i.e.  ensenble  averages  are  equal  to 
the  corresponding  time  averages  of  one  of  the  repre- 
sentative members  of  the  ensemble. 

(3)  The  field  vectors  associated  with  the  source  possess 
fourier  transforms.^ 


The  real  vector  potential  of  the  source  at  a point  in 
space  S (the  bar  under  the  symbol  will  designate  a vector 
quantity)  and  time  t is  A(S,t).  In  the  frequency  domain  the 
fourier  transform  of  A(S,t)  can  be  written: 


A (S , f ) 


r A(S,t)£.i27Tftdt 

V CO 


II-l 


and  from  the  theory  of  fourier  integrals 
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one  has: 


A (S , t)  = 


A(S,f)e 


-i2nf  t 


df 


II-2 


Equations  II-l  and  II-2  define  a fourier  transform  pair  A 
and  A.  Next,  denote  the  time  average  of  a quantity  y(t) 
as  and  define  ^ y ^ as:"* 


LimitJ_  fT 

/y\=  T-> °°2T  l y(t)dt 
V_J -T 


In  most  cases,  y(t)  is  given  by 


II-3 
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y (t)  |t|<  t;] 

0 1 1 1 > T. 


Then  I I- 3 becomes 


Limit 
= T-»  T 


_1_/~T 

2T\  y (t)dt 


II-3a 


The  average  intensity  of  the  radiation  at  a point  S 


is  defined  as 


I(S)  = <^jA(S,t)  |^). 


II-4 


To  find  the  intensity  of  the  disturbance  at  the  detector, 
consider  the  vector  potential  A incident  on  the  beamsplitter- 
compensator.  The  vector  potential  is  divided  into  two  parts 
and  A^  that  travel  to  mirrors  and  M2.  Each  disturbance 
makes  one  reflection  and  one  and  one  half  passes  through 
the  beamsplitter-compensator  before  reaching  the  detector. 

Let  the  complex  reflectance  and  transmittance  of  the  beam- 
splitter-compensator  be  r (f ) and  n(f),  and  let  the  difference 
in  time  of  arrival  of  the  two  disturbances,  resulting  from 
a difference  in  path,  be  t.  Then  the  magnitude  of  A^  and  A2 
are  given  by 

A1  (S,t)  = rr)A  (S,  t) 

A2(S,t)  = rr)A  (S , t+x ) 

The  resultant  amplitude  at  the  detector  is  then 


Ares(So'fc)  = rnA(S0,t)  + rnA(S0,t+x)  11-5 

where  SQ  is  the  position  of  the  detector,  and  the  intensity 


at  the  detector,  I(Sq),  is  given  by 


<so>  • <jAreSlso't)l2^> 


11-6 


r 


7 


from  equation  II-4.  Using  equation  II-5,  and  property  (1) 


W 

V - 

ri 


of  the  source,  equation  11^6  becomes 

I (SQ)  = 2 | rn | 2 


|"<fA(S0,t)  |2)  + <A(S0, 
mated,  one 
<^A(S0,t)  |^>  = ^A(S,t)  |^=  I 


t) "A (S_ , t + t] 


Since  the  source  is  collimated,  one  can  write 


where  I denotes  the  intensity  of  the  incoming  radiation. 
Dropping  the  SQ,  equation  II-7  becomes 

I(Det)  = 2 1 r | 2 * YI  +^(t)  * A(t+xj>>  11-8 

Define^ 

r 11  ( T ) = <A(t)  • A(t+X)^> 

as  the  auto-correlation  function  of  A(t) , and  define^ 

Yh(t)  = ru(T)/I 

as  the  normalized  auto-correlation  function  of  A(t).  Then 
equation  II-8  becomes:  * 

I (Detector)  = 2 | rn | 2 I *^1  + Yh(t)J  II-9 

for  a time  difference  x.  Equation  II-9  is  a form  of  the 
basic  equation  used  in  FTS  and  expresses  the  intensity 
at  the  detector  as  a result  of  a disturbance  with  intensity 
I incident  on  the  beamsplitter-compensator,  and  a time  dif- 
ference x resulting  from  the  radiation  traveling  a path 
difference  X in  the  two  arms  of  the  interferometer.  To 
express  equation  II-9  as  a function  of  path  difference  and 
the  power  spectrum  of  the  source,  consider  equation  II-4. 
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II-7 


J 


Substituting  the  expression  for  A from  equation  II-2  into 


II-4  one  finds: 


Limit  1 


I(S)  = T — 7 00  2T  |A.(S,t)  pdt 


Limit  1 


11-10 


= T — > 00  2T  A (S , f ) df 
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which  is  essentially  Parseval's  theorem.  Define  ' ' 


the  following 

Limit  | A (S , f ) | 2 

G(S,f)  = T — 2T 


11-11 


as  the  power  spectrum  in  the  frequency  domain  f.  Since 


A (S , t)  is  real  it  can  be  shown  that  equation  11-10  can 


be  written: 


11-12 


I(S)  = 2 \ G (S  , f ) df 


and  it  can  also  be  shown  that  the  normalized  auto-corre- 


lation  function  (dropping  the  spatial  dependence)  can 


be  written 


♦ G (f ) cos2irf  tdf 


11-13 


G (f ) df 


Substituting  11-12  and  11-13  into  II-9  gives  (I) is  the 
intensity  at  the  detector 


I(t)  = 2|rn| 


2G(f)  [l+cos2Trf t]df . 


11-14 


Finally,  going  from  the  frequency- time  domains 
to  wavenumber-length  domains  one  makes  the  following 
substitutions : 
x 

t = c c = speed  of  light 


f = co  o = wavenumber (cm-1) 

and  let  2cG(co)  = B(a)be  the  power  spectrum  in  the  wavenumber 
domain,  then  equation  11-14  becomes 


I (x)  = 2 | rn| 


B (o)  [l+cos2irox]  do 


11-15 


Where  x = path  difference  between  the  two  mirrors 

r = reflectance  of  the  beamsplitter-compensator 

n = transmittance  of  the  beamsplitter-compensator 

B(o)  = power  spectrum  of  the  radiation  in  the  wavenumber 
domain 

I (x)  = the  intensity  at  the  detector  for  a path 
difference  x 

Thus  I (x)  consists  of  a constant  term  plus  a fluctuation 
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term.  Let  the  fluctuating  part  of  I(x)  be  dl(x)  then 
equation  11-15  becomes 


dl  (x)  = 2|rn|2  * B (o)cos2Ttoxdo.  11-16 

From  the  orthogonality  properties  of  the  cosine  functions, 
equation  11-16  can  be  solved  for  B(o)  and  one  finds 

B(o)  = 

Equations  11-16  and  11-17  are  the  basic  equations  of  FTS. 

The  quantity  dl (x)  is  called  the  interf erogram,  and  equation 
11-17  shows  that  the  power  spectrum  of  the  source  can  be 
obtained  from  a fourier  cosine  transform  of  the  interfero- 
gram.  It  is  also  clear  from  equation  11-16  that  the  intensity 
at  the  detector,  for  a given  path  difference,  contains 
information  about  the  entire  spectral  region  of  B(o). 

C.  Recovery  of  the  Spectrum 

Implementation  of  the  FTS  technique  requires  the  inter- 
ferogram  dl (x)  to  be  recorded  as  a function  of  path  difference 
x.  In  practice,  x is  continuously  varied  by  moving  one  of 
the  mirrors  at  velocity  v while  the  fluctuating  component 
of  the  signal  from  the  detector  is  recorded.  Two  difficulties 
immediately  arise  in  the  actual  measurement. 


1)  In  general,  equation  11-17  must  be  solved  by  a 
digital  computer  so  dl (x)  must  be  digitized  and 
recorded  as  a discreet  set  of  points  dl(Xj)  as 
opposed  to  an  analytic  function. 

2)  The  limits  on  equation  11-17  imply  that  one  of  the 
mirrors  must  be  moved  an  infinite  distance.  Since 
this  is  not  possible,  the  mirror  is  moved  only  over 
some  finite  distance  L. 


Consider  the  effects  of  these  departures  from  the 
ideal  theoretical  equations  11-16  and  11-17.  To  arrive 
at  the  effect  of  finite  mirror  travel  define  the  function 

1 x < x 
0 x > x 

where  X is  the  maximum  path  difference  attainable.  Equation 
11-17  can  be  written 


a (x)  = 


^Vx> 


Br (a)  = |rn 


a (x)  dl  (x)  cos27roxdx 


where  B^fo)  denotes  the  recovered  spectrum  under  the  conditions 
of  finite  path  difference.  Taking  the  fourier  transform  of 
B r(o)  one  has 

FT[Br(a)]  = a(x)dl(x) 
and  for  the  ideal  spectrum 

FT [B (o) ] = dl (x) 

Then  writing  a(x)  as: 

a (x)  = FT  [a  (o)  ] 
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one  has: 


FT[Br(o)]  = FT [B (a) ] • FT [a (o ) ] . 

Using  the  convolution  theorem®*  one  has  from  the  above; 


Br(o)  = B(o)  * a (a) 


11-18 


where  B(a)  * a (a)  is  the  convolution  of  B(a)  and  a (a) 
and  is  defined  as® 


B (o ) * a(o)  = B (u) a (o-u) du. 


11-19 


To  see  the  effect  of  the  convolution  of  B and  a consider 
the  expression  for  a(o); 


a (a)  = 


sin  2ttqX 

cos27rxadx  = 2t:o 


Thus  one  sees  that  the  recovered  spectrum  Br  is  a smeared 

out  version  of  B due  to  the  convolution  of  B with  the 
2 

sine  function.  The  Rayliegh  criterion  for  resolving 
two  spectral  features  is  that  the  peak  of  one  feature 
falls  at  the  first  zero  of  the  second  feature.  The  first 
zero  of  the  sine  function  occurs  at  2ttcjX  = u,  so  using 
the  Rayliegh  criterion,  the  resolution  6o  is  given  by: 

2n6aX  = n 


6a  = 2' 


11-21 


Thus  the  effect  of  finite  mirror  travel  is  to  determine 
the  resolution  of  the  instrument. 
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The  effects  on  the  recovered  spectrum  due  to  the 
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discreet  nature  of  the  sampled  interf erogram  range  from 
signal-to-noise  degradation^  to  determining  the  spectral 
bandwidth  that  can  be  recovered  from  the  sampled  interfero- 
gram. 2,3,4  It  is  this  latter  effect  that  must  be  carefully 
controlled  in  order  to  avoid  recovering  a false  spectrum. 

The  results  of  sampling  the  interf erogram  will  not  be 
derived  but  will  be  stated.  The  sampling  theorem2,3'4 
states  that  one  must  sample  the  data  at  a rate  equal  to 
the  reciprocal  of  twice  the  highest  frequency  present  in 
the  data  in  order  to  avoid  overlapping  or  aliasing.  In 
equation  11-16,  for  a given  displacement,  the  highest  fre- 
quency present  is  a , where  o is  the  highest  wavenumber 
present  in  the  spectrum  B(o).  Therefore  if  Ax  is  the 
sample  spacing  of  the  interferogram,  the  sampling  theorem 
requires  that: 


1 


Ax  < 2omax  . 11-22 

That  is,  the  sample  spacing  must  be  less  than  or  equal  to 
the  reciprocal  of  twice  the  highest  wavenumber  present  in 
the  spectrum.  If  relation  11-22  is  not  adhered  to,  the 
recovered  spectrum  Br(o)  will  be  aliased  or  folded  back 
on  itself  and  thus  will  not  give  an  accurate  measure  of 
B(o).  In  practice,  one  usually  uses  a low  pass  spectral 
filter  to  insure  that  no  wavenumber  greater  than  omax 


i 
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enters  the  interferometer  for  the  sample  spacing  chosen. 
For  example,  if  one  desires  to  investigate  the  spectrum 
from  a source  in  the  wavenumber  range  0 to  7000cm  \ the 
interferogram  dl£c)  must  be  sampled  every  7.1428  X 10  “’em 
of  patli  difference. 

To  recover  the  spectrum  Br (a)  from  the  interferogram 

dl  (x) , the  procedure  is  thus  to  record  the  interferogram 

dl(x)  at  path  differences  of  x . , where  x.  - x.  , = Ax 

3 3 3~l 

is  given  by  relation  11-22,  over  a path  difference  x 
determined  by  the  desired  resolution  as  given  in  11-21 
and  then  solve  equation  11-17,  which  for  the  discreet  case 
becomes : 

2 

Br  (°n)  = | rn 2 ) = x I dl  (x^)  a (x^  cos2TTOnxin  11-23 

m=0 


where 


xm  = mAx 


« 


and  a (xm)  is  the  truncating  function  described  earlier. 

There  are  other,  more  subtle  steps  that  must  be  taken  before 
Br(o)  can  be  realized,  and  indeed  there  are  subtle  difficulties 
associated  with  the  solution  of  11-23  involving  manipulating 
large  data  files  and  computing  times,  however  the  basic 
procedures  for  recovering  the  power  spectrum  have  been 
outlined.  The  reader  is  refered  to  references  2,3,  and  4 
for  detailed  discussions  of  the  FTS  technique. 
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D.  Advantages  of  FTS 


In  the  introduction  to  this  section  the  multiplex 
nature  of  FTS  was  shown  to  yield  an  advantage  in  either 
short  measurement  time  or  signal-to-noise  over  a sequential 
measurement  technique  such  as  an  ordinary  monochrometer. 

There  are  two  other  distinct  advantages  of  FTS  that  we 
will  briefly  discuss  at  this  time.  As  mentioned  in  the 
previous  section,  the  sample  spacing  of  the  data  points 
in  the  interferogram  is  important  in  determining  the 
maximum  spectral  bandwidth  that  can  be  unambiguously 
recovered  in  FTS.  From  the  example  given,  to  do  mid-infrared 
spectroscopy  (400-4000cm-'1‘)  the  path  difference  (and  position 
of  the  moving  mirror)  must  be  known  very  accurately.  At 
first  this  might  appear  to  be  a great  disadvantage,  and 
indeed  if  the  problem  had  not  been  solved  it  would  be, 
but  the  solution  yields  an  advantage  in  FTS  over  most 
other  spectroscopic  techniques.  To  solve  the  problem  of 
knowing  the  path  difference  very  accurately  the  interfero- 
grams  produced  by  a very  broad  band  source  and  a near 
monochromatic  source  are  used  as  measuring  tools.  It  is 
seen  from  equation  11-16,  by  letting  B(o)be  given  by: 

constant,  o < or 


B (a)  = 


max 


, o > a 


max 
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that  the  interferogram  essentially  becomes  a sine  function 
and  thus  peaks  sharply  at  x=0.  The  broader  the 
spectrum,  the  narrower  the  peak.  In  this  way,  a very  broad- 
band source  can  give  a good  indication  of  zero  path  difference 
for  the  interferometer.  Similarly,  a near  monochromatic 
source  yields  an  interferogram  that  is  essentially  a cosine 
function  of  2nooX,  where  o0  is  the  wavenumber  of  the  source. 
Thus,  a near  monochromatic  source  yields  a waveform  that  can 
be  used  to  measure  *,  provided  0o  is  known.  The  inter- 
ferometer used  in  this  work  actually  consists  of  three 
interferometers  to  utilize  the  interferograms  just  discussed. 
Figure  II-2  shows  the  configuration  for  this  system  of  inter- 
ferometers. The  "white  light"  source  is  a broad-band  incan- 
descant  bulb  that  is  used  as  a reference  point  for  zero  path 
difference,  and  the  HeNe  laser  is  a near  monochromatic 
source  that  supplies  a cosine  wave  by  which  the  path  difference 
can  be  very  accurately  measured.  This  arrangement  yields  an 
absolute  vacuum  wavenumber  determination  to  .01cm  ^ . Normally, 
one  would  be  required  to  use  infrared  standards  and  correct 
for  wavenumber  in  air  to  achieve  this  kind  of  accuracy; 
however,  with  the  FTS  system,  the  accuracy  is  automatically 
achieved. ^ This  advantage  is  sometimes  called  Connes ' 
advantage. 

Another  advantage  of  FTS  is  in  its'  high  throughput 
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or  etendue.  The  energy  available  at  the  detector  is  given 
by: 


I (o)  = T B (a)  E 11-24 

where  T is  the  transmission  of  the  system,  B(o)  is  the 
intensity  of  the  source  at  wavenumber  o,  and  E is  the 
throughput  of  the  system  defined  as; 


E = Aft  11-25 

where  A is  the  usable  cross  sectional  area  of  the  collimated 
beam  and  ft  is  the  solid  angle  subtended  by  the  detector 
collection  optics.  It  can  be  shown2,3,4  that  the  maximum 
allowable  detector  solid  angle  is  governed  by  the  desired 
resolution  and  is  given  by 

2 7T 

ft  = T H-26 

where  R is  the  resolving  power  a /6a.  Using  Equation  11-26, 
equation  11-24  can  be  written 

TB2ttA 
I = R 


or 

IR  = 2ttTBA  11-27 

From  11-27  one  can  see  that  there  is  a trade-off  between 
resolving  power  and  intensity  at  the  detector;  however, 
the  product  IR  does  not  vary.  For  a conventional  spectro- 
meter the  solid  angle  ft  is  given  by2 

wl 

ft  = f2  11-28 
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where  1 is  the  slit  height,  w is  the  slit  width  and  f is 

the  focal  length  of  the  collimator.  If  the  spectrometer 

has  a dispersion  d0/do  the  resolution  is  given  by 

w d£ 

60  = f do 

thus  the  resolving  power  is  given  by 


a d_0  o 

Jo  = R = Do  f w 

so  that  the  slit  width  is 

d6  a “ 

w = do  f R • 11-29 

From  11-24,  11-25,  11-28,  and  11-29  the  product  IR  for  the 

conventional  spectrometer  becomes 

d6  1 

IR(conv)  = TBA  do  of-  Ix-30 

For  a conventional  spectrometer,  a grating  yields  the  highest 

t~h  5 

value  for  the  dispersion  and  in  M*"  order  one  has 

m 


o = f esc  9 


so  that 


de 

do 


= tan  0 


and  for  a grating  angle  of  about  45°  one  has 

= 1 


d0_ 

do 


* * 

Then  equation  11-31  becomes 


IR(conv)  = TBA  f 
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For  a conventional  spectrometer  the  values  (1/f)  are  very 

2 

low  compared  to  2tt  and  typically  are  of  the  order  1/50. 
Thus  FTS  enjoys  a throughput  gain  of  the  order  of  50  over 
a conventional  spectrometer.  This  advantage  is  usually 
called  Jacquinot's  advantage. 

One  other  advantage  of  FTS  is  its  relative  insensi- 
tivity to  stray  light.  Since  each  wavenumber  o is  modu- 
lated at  a different  frequency,  and  only  light  nearly 
parallel  to  the  moving  mirror  is  modulated,  FTS  offers 
a vast  improvement  over  a conventional  spectrometer  in 

4 

insensitivity  to  stray  light.  In  summary,  FTS  offers  at 
least  four  advantages  over  the  conventional  infrared 
spectrometer.  They  are: 

1)  Multiplex  advantage  (Fellgett  advantage) 

2)  High  wavenumber  accuracy  (Connes  advantage) 

I 

3)  High  entendue  or  throughput  (Jacquinot's  advantage) 

4)  Stray  light  insensitivity 

With  fast  detectors  and  high  mirror  velocities,  FTS  also 
has  the  advantage  of  very  rapid  data  acquisition,  which 
along  with  the  above  mentioned  advantages,  make  it  very 
suitable  for  mid-infrared  emission  studies. 
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III.  Description  of  the  Experiment 
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Figure  III-l  shows  the  basic  layout  used  to  study 
the  infrared  emission  from  a small  combustor.  The  major 
components  in  the  experiment  are  A.  the  combustor,  B. 
collection  optics,  C.  interferometer-detector  system,  and 
D.  data  reduction  system. 

A.  The  Combustor 

The  combustor  chosen  for  the  study  was  a flat  flame 

10  _ 
burner  , so  named  for  the  appearance  of  the  flame  it  pro- 
duces. The  flame  has  the  characteristic  that  the  inter- 
action volume  is  large  in  comparison  to  the  boundary  layer, 
so  that  boundary  layer  effect  can  be  ignored.  The  surface 
of  the  burner  had  a diameter  of  approximately  2.5  cm., 
and  could  be  moved  vertically  by  a precision  screw-adjustment. 
The  burner  was  housed  in  a water  cooled,  copper  cylinder  of 
approximately  7.5  cm.  diameter,  equipped  with  ports  for 
observing  the  flame.  The  fuel  used  was  low  purity  methane 
(commercial  grade,  approx.  90%)  which  was  premixed  with  air 
before  being  burned.  Normally,  the  ratio  of  air  to  methane 
was  approximately  7:1  by  volume,  which  yielded  a very  stable, 
flat  flame  that  appeared  to  extend  a few  millimeters  above 
the  burner  surface. 

B.  Collection  Optics 

A requirement  for  the  interferometer  is  that  the  light 
source  be  collimated.  The  collimating  lens  used  was  a 
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100  mm.  focal  length,  f/2  As2s3  " CaF2  separated  doublet 
with  a useful  transmission  range  from  approximately  1250cm  ^ 
to  12500cm- ^ (8pm  to  .8pm).  The  lens  was  free  of  chromatic 
aberration  between  approximately  1800cm  ^ and  2900cm  1 
(5.5pm  to  3.5pm),  and  had  minimum  chromatic  aberration 
throughout  its  useful  range.  A 100mm  focal  length,  f/2, 

NaCl  lens  was  used  briefly  for  measurements  below  1250cm  ^ 
(above  8pm).  As  shown  in  equation  11-26,  the  resolving  power 
of  the  interferometer  is  dependent  on  the  solid  angle  sub- 
tended by  the  detector  collection  optics.  Since  the  flame 
was  in  no  way  a point  source,  and  it  was  impossible  to  place 
an  aperture  in  the  flame,  it  was  necessary  to  limit  the 
effective  source  size  by  some  other  means.  This  was  achieved 
by  placing  approximately  a 1cm  aperture  at  the  entrance 
to  the  interferometer,  45cm  from  the  collimating  lens.  This 
alone,  however,  was  not  sufficient  to  achieve  .06  cm  ^ reso- 
lution, so  an  adjustment  to  the  interferometer  was  necessary 
which  will  be  discussed  in  the  next  part  of  this  section  of 
the  report. 
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C.  The  Interferometer-Detector  System 

The  interferometer  used  was  essentially  an  EOCOM  model 
FMS  7101  with  the  magnetic  tape  system  disconnected.  The 
instrument  has  greater  than  .06cm-^  resolution  capability 
over  the  wavenumber  range  400cm-^  to  7000cm  ^ and  is  capable 
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of  mirror  velocities  up  to  4cm/sec.  Thus,  rapid  scanning 
of  a broad  spectral  retion  is  possible,  with  all  of  the 
advantages  of  the  FTS  technique  discussed  in  section  III. 

The  beamsplitter-compensator  normally  used  was  Fe20^  coated 
CaF2  with  a spectral  range  from  1250cm-1  to  about  7000cm-1 
(8um  to  1.4ym);  however,  a Ge  coated  KBr  beamsplitter- 
compensator  was  briefly  used  for  measurements  below  1250cm-1 
(above  8ym) . Two  detectors  were  used  during  the  course  of 
the  investigation.  Initially,  a pyroelectric  detector 
(triglycine  sulfate,  TGS)  was  used  with  the  CaF2  beamsplitter; 

O 

however,  the  TGS  detector's  poor  detectivity  (D*=:6X10  cm 
. 5 

Hz"  /watt)  and  slow  time  constant  made  spatially  resolved 
measurements  practically  impossible.  Late  in  the  investi- 
gation, a thermoelectrically  cooled  PbSe  detector  was  made 
available  with  a much  higher  detectivity  (D*  = 2X10^)  , a 
much  faster  time  constant,  but  with  a somewhat  more  limited 
useful  spectral  range.  This  allowed  spatially  resolved  data 
to  be  acquired  in  the  wavenumber  range  1700cm-1  to  7000cm  1 
(5.9ym  to  1.4ym).  The  final  beamsplitter-detector  combinations 
chosen  were  the  CaF2  - PbSe  for  the  region  1700cm-1  to  7000cm-1 
and  the  KBr  - TGS  for  the  region  from  400cm-1  to  1700cm  ^ . 

Very  little  data  was  taken  in  the  latter  region  due  to  the 
poor  performance  of  the  TGS  detector.  The  detectors  were 
placed  30°  off  axis  at  the  focal  point  of  a 5cm  diameter, 

10  cm  focal  length  parabola,  to  effectively  increase  their 
collection  area. 
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As  mentioned  in  part  B,  it  was  necessary  to  further 
reduce  the  solid  angle  subtended  by  the  collection  optics 
from  the  limit  set  by  the  aperture  placed  at  the  interfero- 
meter entrance.  To  achieve  .06cm  ^ resolution  at  7000cm  ^ 
the  solid  angle  must  be  less  than  5.39X10  ^ sterradians, 
which  corresponds  to  an  angle  of  deviation  for  rays  entering 
the  interferometer  of  less  than  .237°  or  4.14  mrad.  This 
was  difficult  to  achieve  without  sacrificing  a large  portion 
of  signal,  and  in  view  of  the  fact  that  most  of  the  emission 
from  the  flame  occurred  in  the  vicinity  of  2000cm  ^ , a 
compromise  was  reached.  The  solid  angle  was  limited  to  about 

_*4 

1.8  X 10  sterradians,  which  was  adequate  to  yield  better 

than  .06cm  1 resolution,  and  line  position  accuracy  to  better 

than  . 03cm-^  at  2000cm  \ This  was  achieved  by  limiting  the 

effective  collection  area  of  the  detectors  from  the  =20cm,t 

2 

of  the  parabola  to  about  0.2cm  . This  limiting  of  off-axis 
rays  from  entering  the  interferometer  and  being  detected 
resulted  in  spatial  resolution  in  the  flame  of  approximately 
1.5mm  perpendicular  to  the  optical  axis. 

D.  Data  Reduction  System 

The  data  reduction  system  consisted  of  a Nicolet  1180 
signal  averager/data  processor  with  a CRT  display  for  visual 
inspection  and  manipulation  of  data  and  a Zeta  model  160 
digital  plotter  for  producing  hard  copy  data  plots.  In 
addition,  a Diablo  model  44  dual  disk  random  access  storage 
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device  was  used  for  storing  data  and  machine  code  instruc- 
tions. The  system  was  used  in  two  basic  modes.  For  moderate 
to  low  resolution  6 ajnLcm-1)  , data  was  acquired  into  the  1180, 
processed  (fourier  transformed) , and  the  resulting  spectrum 
stored  on  disk.  This  mode  allowed  data  points  to  be  taken 
at  a rate  of  50KHz,  so  that  a 1cm-1  resolution  interferogram 
(16384  data  points)  could  be  acquired  in  approximately  0.33 
sec.  and  could  be  transformed  in  another  22  sec.  (32K 
transform  points)  . Low  resolution,  4cm--*-  scans  could  be  taken 
in  less  than  0.09  sec.  and  could  be  transformed  in  5 sec.  It 
was  therefore  possible  to  obtain  a spectrum  of  the  burner 
over  the  region  1700cm--*-  to  7000cm  1 (5.9pm  to  1.4ym)  with 
4cm_1  resolution  is  5 sec.  The  second  mode  of  data  collection 
was  necessary  for  higher  resolution  spectra  and  consisted 
of  acquiring  data  directly  onto  disk,  processed,  and  the 
resulting  spectrum  stored  back  on  disk.  This  mode  dictated 
a somewhat  slower  data  collection  rate  of  20KHz,  so  that  a 
.5cm-1  resolution  scan  (32768  data  points)  could  be  acquired 
in  1.64  sec.  and  could  be  transformed  (64K  transform  points) 
in  about  3 minutes.  A .06cm-1  scan  (262000  data  points) 
required  13  sec.  to  complete,  and  required  about  28  minutes 
to  process  (524K  transform  points) . In  comparison,  the  time 
required  to  duplicate  the  . 06cm'1  resolution  scan  of  the  same 
spectral  region  with  a conventional  monochrometer  (with  the 


I 
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same  signal  to  noise  ratio)  would  have  required  322  hours 
(assuming  the  throughput  of  the  two  spectrometers  is  the 
same)  and  would  still  not  yield  the  stray  light  advantage, 
or  Connes'  advantage  (see  section  II).  An  important  feature 
of  the  1180  data  reduction  system  was  its  capability  for 
co-adding  scans  from  the  interferometer.  In  this  way,  the 
signal  to  noise  ratio  of  the  final  spectrum  could  be  improved 
by  as  much  as  a factor  5000  over  the  signal  to  noise  ratio  of 
a single  scan.  This,  of  course,  required  longer  measurement 
times  so  usually  some  compromise  was  reached  concerning 
acceptable  signal  to  noise  requirements.  The  data  reduction 
system  also  allowed  for  data  file  manipulation  (subtracting 
spectra,  wavenumber  to  wavelength  conversion,  scaling,  etc.) 
which  greatly  enhanced  the  capabilities  of  the  interferometer. 


IV.  Results 


Emission  spectra  were  obtained  of  the  flat  flame  burner, 
using  a pre-mixed  CH4~air  mixture,  in  the  spectral  region 
1700cm-1  to  7000cm  1 (5.9ym  to  1.4ym)  with  approximately 
a 1.5  meter  air  path  between  burner  and  detector.  Spectra 
were  obtained  under  a variety  of  experimental  conditions 
including  variation  of  the  fuel-air  mix  ratio,  moderate  and 
high “resolution,  and  collecting  data  from  different  regions 
of  the  flame.  Some  preliminary  results  were  also  obtained 
in  the  spectral  region  from  400cm  1 to  1700cm  1 (25ym  to 
5.9ym);  however,  data  was  only  obtained  under  one  set  of  the 
aforementioned  conditions.  Only  relative  intensities  were 
obtained,  therefore  only  a wavenumber  scale  appears  on  the 
data  plots,  the  vertical  scale  being  arbitrary  units.  The 
data  was  not  corrected  for  the  spectral  response  of  the  system, 
due  to  the  late  acquisition  of  the  PbSe  detector;  however,  a 
brief  discussion  of  the  general  behavior  of  the  system  function 
will  be  given  later.  The  data  plots  will  be  referred  to  in 
the  text  by  figure  number,  i.e.  figure  IV-1,  however,  many 
of  the  plots  occupy  more  than  one  report  page  so  these  figures 
will  be  labeled  IV-1.1,  IV-1. 2,  etc.  to  indicate  multi-page 
figures . 


A.  Survey  Spectra  - Since  the  flame  was  observed  through 
air,  it  was  necessary  to  realize  the  effects  of  infrared 
absorption  by  the  atmosphere  in  the  spectral  range  of  interest. 
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Figure  IV-1  is  the  spectrum  of  a 1500 °K  "glowbar"  measured 
through  1.5  meters  of  air  with  approximately  1cm""1  reso- 
lution. Data  was  collected  up  to  7000cm  1;  however,  no 
absorption  was  observed  above  5000cm  1 so  the  plot  was  ter- 
minated at  this  point.  The  glowbar  spectrum  is  roughly  that 
of  a black-body  radiator,  therefore  it  produces  a continuous 
spectrum  that  can  be  approximated  by  the  Planck  equation.1® 

The  collection  optics,  detector,  and  beamsplitter  were 
as  discussed  in  section  III.  The  number  of  scans  on  the 
plot  indicates  the  number  of  single  scan  interferograms 
that  were  co-added  before  being  processed.  It  can  be  seen 
that  there  is  strong  absorption  in  three  regions  of  the 
spectrum;  a)  water  vapor  absorption  in  the  vicinity  of 
1700cm-1  (6.3ym  band),  b)  C02  absorption  in  the  vicinity  of 
2350cm  1 (4.3ym  band),  and  c)  water  vapor  absorption  ( 2.7pm 
band)  and  C02  absorption  (2.7ym  band)  in  the  vicinity  of 

3700cm-1.11  There  is  also  some  weak  absorption  in  the  vicinity 

-1  . 11 

of  3100cm  which  is  also  due  to  water  vapor.  By  comparing 

figure  IV-1  with  a theoretical  black-body  function,  the 

approximate  behavior  of  the  system  response  function  can 

be  obtained.  The  response  function  obtained  in  this  manner 

is  approximate  because,  a)  the  glowbar  does  not  yield  a true 

black-body  spectrum,  and  b)  the  presence  of  atmospheric 

absorbers  distorts  the  measured  spectrum.  The  general 
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behavior  of  the  system  response  function,  however,  can  be 
obtained  in  this  manner.  Figure  IV-2  shows  the  results  of 
this  calculation.  The  curve  was  obtained  by  dividing  points 
on  the  measured  plot  by  the  predicted  irradiance  from  a 
black-body,  normalizing  the  data  to  the  maximum  value, 
and  interpolating  the  data  with  a cubic  spline.  The  sharp 
drop  in  sensitivity  below  1950cm_^  is  due  to  the  cutoff  of 
the  PbSe  detector  and  the  broad  dip  in  the  curve  between 
2600cm  1 and  3000cm  ^ is  due  to  the  transmission  characteristics 
of  the  As2S.j  lens  component.  The  gradual  decrease  in  sensi- 
tivity above  2200cm-^  is  also  primarily  the  result  of  the 
decreased  responsivity  of  the  PbSe  detector. 

Figure  IV-3  is  the  spectrum  obtained  of  the  flat  flame 
burner  with  approximately  0.25  cm  ^ resolution.  As  in 
figure  IV-1,  data  was  obtained  from  1700cm  1 to  7000cm  "S 
however,  no  interesting  features  were  observed  above  5000cm  ^ 
so  the  plot  was  terminated  at  this  point.  The  data  was 
collected  from  a region  of  the  flame  approximately  10mm  above 
the  burner  base.  The  actual  shape  of  the  region  over  which 
data  was  collected  is  difficult  to  calculate  precisely; 
however,  from  simple  ray  tracing  one  finds  that  the  volume 
sampled  is  roughly  a cone  11mm  high  and  1.5mm  at  the  base, 
oriented  along  the  optical  axis  of  the  system.  Since  the 
flame  is  very  uniform  along  the  optical  axis,  the  important 
parameter  is  the  height  above  the  burner  base  the  data  was 
taken.  The  strongest  features  in  the  spectrum  are  a)  pre- 
dominantly H2O  emission  between  1700cm  ^ and  2050cm  ^ , 
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b)  CO^  emission  from  2050cm  * to  2400cm  1 c)  H20  and  CC>2 

from  3200cm-'''  to  4100cm-1.  12  Figures  IV-1  and  IV-2  should 

be  kept  in  mind  while  considering  figure  IV-3.  The  H20 

emission  between  1700cm-1  and  2050cm-1  occurs  in  a region  of 

relatively  strong  atmosphere  absorption  and  relatively  poor 

detector  response,  therefore  the  intensity  of  the  band  is 

considerably  greater  than  it  appears  on  the  plot.  Similar 

arguments  apply  for  the  H20  and  CC>2  bands  from  3200cm  1 to 

4100cm-1.  Atmospheric  absorption  is  also  present  in  the 

region  from  2300cm-1  to  2380cm  ^ in  the  region  of  CC>2  emission 

however,  in  the  regions  from  2050cm  1 to  2300cm  1 and  from 

2380cm-1  to  2400cm  1 there  is  little  absorption.  The  lower 

frequency  H20  emission  results  from  transitions  involving 

a change  of  one  quantum  of  the  V2  normal  mode  frequency, 

i.e.  Av2=1.11  The  CC>2  emission  in  the  vicinity  of  2300cm-1 

results  from  transitions  involving  a change  of  one  quantum 

11  12 

of  asymetric  stretch,  i.e.  Av2=l.  ' The  density  of  lines 
in  this  region  is  so  great  that  the  line  spacing  is  less  than 
the  doppler  width  of  the  lines,  therefore  the  feature 
appears  to  be  almost  a continium  with  fine  structure  super- 
imposed on  it.  Two  band  heads  can  be  seen  at  2384  cm  1 
and  2397cm-1  resulting  from  the  011!  Or^O  and  00°1  -»■  00°0 
transitions.  The  H20  and  C02  emission  in  the  high  frequency 
region  results  from  a number  of  transitions;  however,  the 
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stronger  one  for  H^O  are  Av^=l,  Av3=l,  and  AV2=2  and  for 
CO2,  the  combination  bands  Av^=l  + Av^  = 1 and  Av2  = 2 + 


B.  Emission  from  Different  Zones  - Figures  IV-4  through 

IV- 8 show  the  spectra  obtained  at  various  heights  above 

the  burner  base.  The  emitting  region  has  a height  of  about 

1.5mm,  except  for  the  region  at  the  burner  base,  where  the 

base  acts  as  an  additional  stop  in  the  collection  optics, 

reducing  the  height  of  the  region  to  about  .75mm.  The  CH  : 

4 

Air  ratio  was  1:7  by  volume  in  all  of  the  data.  Only  the 
spectral  region  from  1700cm  1 to  2400cm  1 is  plotted,  as  the 
remaining  regions  showed  little  change.  The  resolution  in 
all  the  plots  is  1cm  1 . While  the  vertical  scale  is  arbi- 
trary, comparisons  can  be  made  between  plots  in  this  series 
as  factors  such  as  system  gain,  number  scans,  etc.  was  constant. 
Due  to  the  smaller  emitting  volume  at  the  burner  base,  the 
vertical  scale  in  figure  IV-4  is  expanded  by  2 over  the  re- 
maining plots.  An  interesting  feature  in  the  plots  is  the 
appearance  of  lines  from  the  CO  fundamental  in  the  data  near 
the  base  of  the  burner.  The  structure  extends  from  about  2000 

cm-1  to  2250cm-1  and  results  primarily  from  the  1 -*■  0 transi- 
14 

tion.  The  CO  emission  drops  off  very  quickly  above  the  base 
indicating  combustion  is  going  to  completion  very  rapidly. 

This  is  consistent  with  mass  spectrometer  data  from  Battelle 
Laboratories  showing  CO  concentration  going  to  zero  approxi- 
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inately  5mm  above  the  burner  base  for  a 1/10  atmosphere 

15  . 

flame.  At  atmosphere  the  combustion  is  expected  to 

occur  more  rapidly,  thus  the  CO  concentration  should  drop 

off  even  faster  in  our  case.  No  CO  overtone  spectrum  was 

observed  in  any  of  the  data. 

C.  Variation  of  Mix  Ratio  - Figures  IV  - 9 through  IV  - 12 
show  spectra  oblained  with  rich  CH^  - Air  mixtures  ranging 
from  1:3.4  to  1:5.8  by  volume.  All  data  was  taken  at  a 
height  above  the  burner  base  of  10mm.  As  in  part  B) , in- 
tensity comparisons  among  the  plots  can  be  made.  The  reso- 
lution was  1cm-1.  Again,  CO  emission  is  observed  between 
2000cm-1  ana  2050cm-1,  for  mix  ratios  less  than  approximately 
5:1.  In  this  data,  2 -►  1 transitions  can  be  seen  along  with 
the  1 -*■  0.  No  methane  bands  could  be  identified  in  any  of 
the  data,  nor  could  any  CO  overtone  bands  be  seen. 

D.  High  Resolution  Spectra  - Figures  IV  - 13,  14  and  15 
show  spectra  taken  with  .06cm-1  resolution  in  the  region 
1700cm-1  to  2400cm-1  under  two  sets  of  experimental  para- 
meters. Figures  IV  - 13  and  IV  - 14  were  obtained  from  a 
region  of  the  flame  just  above  the  burner  base  with  a "normal" 
1:7  CH^:Air  mix,  while  figure  IV  - 15  was  obtained  at  a height 
of  6mm  above  the  burner  base  with  a rich  1:3.9  CH^Air  mix. 

Figure  IV  - 14  is  a detail  plot  of  the  rotational  structure 

° ° . 
near  the  00  1+  00  0 band  head  of  C02-  By  comparing  the  spec- 
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tra  with  those  of  Pyler,  et  al.  it  is  found  that  J"  values 
up  to  102  are  observed.  In  both  figures  IV  - 13  and  IV  - 15 
it  is  clear  that  the  C02  line  density  in  the  region  from 
roughly  2100cm-1  to  2300cm-1  is  too  great  to  allow  individual 
lines  to  be  resolved  at  this  temperature  and  pressure.  Emis- 
sion from  the  CO  fundamental  band  can  be  seen  in  both  figures. 

E.  Data  Below  1700cm-1  - Some  preliminary  results  were  ob- 
tained in  the  region  from  about  400cm  1 to  1700cm  1/  with 
roughly  1cm  1 resolution,  using  the  TGS  detector  and  the 

KBr  beamsplitter.  A sample  of  the  spectra  is  shown  in  figure 
IV  - 16.  As  discussed  earlier,  the  poor  performance  of  the 
TGS  made  data  aquisition  in  the  region  very  difficult,  thus 
only  one  set  of  parameters  was  used.  The  CH^:Air  mix  was 
1:7  and  the  region  of  the  flame  sampled  was  6mm  above  the 
burner  base.  The  additional  features  below  1700cm  1 are 
premarily  the  15um  band  of  C02  (in  the  vicinity  of  670cm  1) 
and  more  of  the  6.3um  band  of  H20  (in  the  vicinity  of  1600cm  1) . 

F.  Rotational  Temperature  Determination  - The  rotational 
temperature  of  CO  in  a rich  1:3.9  CH^:Air  mixture  was  cal- 
culated from  the  relative  intensities  of  the  rotational  lines 
of  the  fundamental  from  the  data  in  figure  IV  - 15.  The 

intensity  of  a rotational  emission  line  from  a given  vibra- 

. 14 

tional  level  can  be  written  as 
2Nv4 

Ijf  ji(  = -Q^~  Sj„Exp[-F(J')hc/kT]  IV  - 1 
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where  N is  a constant  for  the  given  vibrational  level,  v is 

the  wavenumber  of  the  line,  Qr  is  the  rotational  partition 

function,  F(J')  is  the  rotational  term  value  of  the  upper 

level,  h is  Planck's  constant,  c is  the  speed  of  light,  k is 

Boltzmann's  constant,  and  T is  the  rotational  temperature. 

S is  the  part  of  the  square  of  the  transition  moment  that 
J 

depends  on  J and  is  given  by: 

S j„  = J"  + 1 for  AJ  = +1 
SJM  = J"  for  AJ  = -1. 

The  procedure  used  to  determine  N and  T was  to  fit  the  mea- 
sured relative  intensities  of  a series  of  rotational  lines, 
from  a given  vibrational  level,  to  equation  IV  - 1 using  a 
non-linear  least  squares  adjustment  due  to  Marquardt . ^ 
Therefore,  each  vibrational  band  of  the  fundamental  that 
could  be  observed  yielded  a rotational  temperature.  Due  to 
interference  from  C02  emission,  only  P branch  lines  were 
used  which  ranged  in  wavenumber  from  1972cm  ^ to  2056cm 
Since  the  wavenumber  range  was  limited,  no  correction  was 
made  for  system  sensitivity.  Line  identification  was  made 

by  comparing  measured  lines  with  lines  calculated  from  the 

17  ... 

constants  of  Mantz,  et  al.  , and  lines  identified  as  con- 
sisting of  unresolved,  multiple  lines  were  not  used  in  the 
calculation.  By  comparing  the  fuel  rich  spectrum  with  a 
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"normal"  mix  spectrum,  CO  lines  which  overlapped  1^0  lines 
were  also  excluded  from  the  calculation.  Thirty  seven 
lines  remained  in  the  P branch  which  were  judged  sufficiently 
interference  free  to  be  used  in  the  calculation.  Figure  IV-17 
is  a sample  of  the  line  identification  procedure  used.  Lines 
were  coarsely  identified  on  the  plot  by  computer,  then  were 
checked  for  interference  from  higher  vibrational  levels 
using  tables  constructed  from  the  constants  of  ref.  17.  Due 
to  the  size  of  the  source,  line  positions  were  only  measured 
to  .03cm-l  accuracy.  Resolution  was  less  than  the  half  widths 
of  the  lines. 

The  temperature  calculated  from  the  three  lower  vibra- 
tional levels  observed  are  given  below. 


o 

II 

> 

7 

T = 2199 

+ 

160 

°K 

V"=l 

7 

T = 1936 

+ 

53 

°K 

V"=2 

7 

T = 1958 

+ 

86 

°K 

The  errors  shown  are  the  estimates  of  the  standard  deviations 
from  the  least  squares  routine.  Since  the  data  was  not  cor- 
rected for  self  absorption,  the  temperature  determined  from 
the  1 ->•  0 vibrational  band  is  expected  to  be  in  error  by  the 
largest  amount,  probably  by  more  than  the  160 °K  estimated. 

The  weighted  average  of  the  three  values  was  calculated  using 
the  estimated  standard  deviations  which  yielded; 

T (weighted  average)  = 1963 °K 
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The  estimated  error  in  the  temperature  was  taken  to  be  the 
average  of  the  upper  two  levels'  standard  deviations  since 
the  error  in  the  V"=0  temperature  is  probably  largely  due 
to  self  absorption  and  not  statistical  fluctuations.  The 
constants  N(0),  N(l),  and  N(2)  were  determined  to  be; 

N ( 0 ) = 5.00E  - 13  + 2E  - 16 

N ( 1 ) = 2.35E  - 13  + 6E  - 17 

N ( 2 ) = 9.60E  - 14  + 6E  - 17 


where  the  units  are  arbitrary.  Using  the  average  rotational 
temperature  and  the  calculated  constants,  the  resultant  line 
intensities  for  the  three  bands  were  calculated  and  are  listed 
in  table  IV- 1.  Also  listed  are  the  calculated  wavenumbers, 
the  measured  intensities,  the  percent  difference  between  cal- 
culated and  measured  intensities,  and  the  line  designations. 

As  mentioned  earlier,  all  lines  are  P branch  lines.  As- 
suming the  rotational  temperature  is  a reasonable  measure  of 
the  gas  kinetic  temperature,  the  average  temperature  agrees 
well  with  measured  values  from  Battell.^ 

It  should  be  noted  that  the  values  of  N(0),  N(l),  and 
N(2)  should  yield  information  about  the  relative  vibrational 
level  populations.  Since  ST„  is  the  J - dependent  part  of 

the  square  of  the  transition  moment,  the  constants  N(v")  can 

, . 14 

be  written. 


N(v") 


W n ( v ' ) 


IV-2 
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Table  IV-1 


Wavenumber 
( cm"  ) 

Designation 

(v",J") 

Intensity  Meas. 
(arb.  units) 

Intensity  Calc, 
(arb.  units) 

Percent  Diff. 

1971.98 

(2,28) 

4 . 210E-02 

4 . 062E-02 

3.51 

1978.31 

(1,32) 

9. 249E-02 

8 . 305E-02 

10.21 

1978.93 

(0,37) 

1.459E-01 

1. 270E-01 

12.92 

1983.13 

(1,31) 

9 . 463E-02 

8 . 848E-02 

6.50 

1983.94 

(0,36) 

1 . 531E-01 

1. 380E-01 

9.88 

1985.88 

(2,25) 

4 . 202E-02 

4 . 617E-02 

9.88 

1987.92 

(1,30) 

9 . 212E-02 

9 . 391E-02 

1.94 

1988.91 

(0,35) 

1. 529E-01 

1.493E-01 

2.34 

1990.45 

(2,24) 

4 . 459E-02 

4 . 777E-02 

7.13 

1993.86 

(0,34) 

1. 733E-01 

1. 610E-01 

7.09 

1994.99 

(2,23) 

5. 344E-02 

4 . 920E-02 

7 . 94 

1997.42 

(1,28) 

9. 566E-02 

1. 045E-02 

9 .28 

1998.78 

(0,33) 

1 . 779E-01 

1. 730E-01 

2.77 

1999.50 

(2,22) 

6 . 233E-02 

5. 043E-02 

19.09 

2002.11 

(1,27) 

1. 065E-01 

1. 096E-01 

2.93 

2003.67 

(0,32) 

2. 035E-01 

1. 851E-01 

9.03 

2003.98 

(2,21) 

5 . 282E-02 

5 . 145E-02 

2.60 

2006.78 

(1,26) 

1. 145E-01 

1. 145E-01 

.02 

2011.42 

(1,25) 

1 . 244E-01 

1 . 190E-01 

4 . 32 

2013.35 

(0,30) 

2. 200E-01 

2. 096E-01 

4.71 

2016.03 

(1,24) 

1 . 220E-01 

1. 232E-01 

1.00 

2017.21 

(2,18) 

4. 800E-02 

5.292E-02 

10.26 

2018.15 

(0,29) 

2 . 263E-01 

2 . 218E-01 

2.00 

2021.56 

(2,17) 

4 . 611E-02 

5. 282E-02 

14.55 

2022.91 

(0,28) 

2 . 378E-01 

2. 337E-01 

1.74 

2025.14 

(1,22) 

1. 299E-01 

1 . 302E-01 

.25 

2027.65 

(0,27) 

2. 376E-01 

2 . 4 52E-01 

3.20 

2029.66 

(1,21) 

1. 310E-01 

1. 329E-01 

1.46 

2032.35 

(0,26) 

2. 514E-01 

2 . 562E-01 

1.91 

2034.41 

(2,14) 

4 . 7 31E-02 

5. 047E-02 

6.68 

2037.03 

(0,25) 

2 . 729E-01 

2. 666E-01 

2.32 

2041.67 

(0,24) 

2 . 731E-01 

2. 761E-01 

1.10 

2043.00 

(1,18) 

1. 401E-01 

1. 369E-01 

2.28 

2046.28 

(0,23) 

2. 763E-01 

2. 847E-01 

3.03 

2046.98 

(2,11) 

5 . 170E-02 

4 . 486E-02 

13.22 

2050.85 

(0,22) 

2. 854E-01 

2 . 921E-01 

2.36 

2056.05 

(1,15) 

1 . 387E-01 

1 . 337E-01 

3.64 
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where  W is  a constant,  n(v')  is  the  number  density  of  mole- 


• , , \r  * y ** 

cules  in  the  v'  vibrational  level,  and  R ' is  the  transi- 


tion moment.  For  a linear  dipole  moment  function,  and  Av=l, 

1 8 

it  can  be  shown  from  the  calculations  of  Heaps  and  Herzberg 
[ v 1 v11!^ 

that  |R  ' | can  be  approximated  by: 

i 2 


R 


v 1 , v 


»l  2 


W, 


R 


1, 


T • 


0^0 
V ' , v" / 


if  2v 1 <<  — tt  - 1 . In  equation  IV-3  v , 

WeXc  v , v 


IV-3 

is  the  wave- 


number  for  the  J"=J'=0  transition.  Using  equation  IV-2  and 
IV-3,  the  relative  number  densities  found  are: 


n(2)/n(l)  = .229 
n(3)/n(l)  = . 0609 

Assuming  the  vibrational  levels  follow  a Boltzman  distribu- 
tion, the  previous  ratios  would  yield  a vibrational  tempera- 
ture of; 


T (vibrational)  = 2115  + 50°K. 
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V.  Conclusions 


FTS  can  yield  rapid,  high  resolution  infrared  emission 
spectra  of  combustion  plasmas  over  a broad  spectral  region 
with  high  sensitivity,  high  wavenumber  accuracy,  and  good 
stray  light  rejection.  If  high  resolution  is  not  needed, 
very  rapid  scanning  and  processing  can  be  accomplished. 
Manipulation  of  stored  spectra  (background  subtraction, 
wavenumber  to  wavelength  conversion,  etc.)  can  be  performed 
using  the  same  computing  equipment  needed  for  fourier  trans- 
forming the  interferograms . Once  a particular  optical  sys- 
tem and  detector  have  been  calibrated  with  a standard  black- 
body,  absolute  intensity  measurements  become  routine. 

Species  identification,  remote  temperature  determination,  and 
species  concentration  determination  (if  absolute  intensities 
are  known)  is  possible  if  infrared-active  species  are  pre- 
sent in  the  plasma.  The  technique  overcomes  many  of  the  low 

signal  level  problems  normally  associated  with  infrared 

« 

spectroscopy  and  thus  allows  the  abundant  information  that 
exist  in  the  infrared  region  to  be  more  readily  obtained. 
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